Size  Scale  and  Defect  Engineered  Nanostructures  for  Optimal  Strength  and  Toughness 


Final  Report 

PI:  W.W.  Gerberich,  University  of  Minnesota 


I.  Introduction 

The  goal  of  this  research  project  was  to  develop  realistic  combinations  of 
ceramics  and/or  semiconductors  that  simultaneously  achieve  high  hardness  (>40  GPa) 
and  toughness  (>400  MPa*m*^^).  While  this  is  probably  unrealistic  for  bulk  material 
production  due  to  inherent  defects  associated  with  processing,  there  is  a  very  real 
opportunity  for  thin  coatings  in  environments  involving  severe  wear  and/or  dynamic 
penetration.  This  includes  a  wide  range  of  industrial  applications,  including  automotive, 
aircraft,  electronics,  manufacturing,  and  biomedical.  To  achieve  this  goal,  the  proposal 
outlined  three  primary  objectives:  (1)  the  development  of  physically  based,  dislocation 
models  to  understand  the  deformation  of  brittle  materials,  (2)  the  synthesis  of  model 
ceramic  nanocomposites  that  demonstrate  high  hardness  and  toughness,  and  (3)  the 
detailed  understanding  of  the  arrangements  and  types  of  dislocation  structures  in  small 
volumes.  This  involved  the  uniaxial  compression  of  Si  nanovolumes  (spheres  and 
towers)  using  a  combination  of  TEM  in  situ  indentation  and  molecular  dynamics 
simulations  for  objective  (1),  the  deposition  of  Si-SiC  core-shell  nanotowers  for 
objective  (2),  and  the  HR  TEM  analysis  of  deformed  Si  nanovolumes  for  objective  (3). 
This  report  summarizes  the  substantial  progress  made  on  the  first  two  of  these  objectives 
and  the  work  that  is  currently  underway  to  address  objective  (3). 


II.  Uniaxial  compression  of  Si  nanospheres  and  towers 

One  of  the  major  components  of  this  research  project  was  the  deformation  of 
silicon  nanospheres  and  nanotowers  through  TEM  in  situ  nanoindentation  experiments. 
These  consisted  of  a  series  of  uniaxial  compressions  of  hypersonic  plasma  particle 
deposition  (HPPD)  [1,  2]  synthesized  single  crystal  Si  nanospheres  (J=63-349  nm)  on 
AI2O3  substrates  and  VES  grown  Si(lll)  nanotowers  (J=23 1-415  nm).  Figure  la,b  gives 
an  example  of  a  173  nm  diameter  Si  nanosphere  before  and  after  compression.  In  the 
nanospheres,  the  oxide  shell  prevents  the  release  of  dislocations  and  leads  to  the 
formation  of  a  back  stress  due  to  dislocation  pile-up.  The  applied  shear  stress,  r,  needed 
to  nucleate  additional  dislocations  is  thus  increased  and  can  be  described  by  an  Eshelby 
dislocation  pile-up  mechanism  [3]  as 
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where  fj.  is  the  shear  modulus,  b  is  the  Burgers  vector,  N  is  the  number  of  dislocations,  v 
is  Poisson’s  ratio,  and  d  is  the  sphere  diameter.  Figure  Ic  shows  the  average  flow  stress 
as  a  function  of  diameter  for  nanospheres  found  using  a  geometric  contact  area.  In 
addition,  ex  situ  indentations  of  Si  nanospheres  are  included  from  Mook  et  al.  [4],  where 
the  flow  stress  was  calculated  based  on  an  average  between  a  geometric  contact  and 
cylindrical  contact  [4].  Assuming  =  It  and  using  =  60.5  GPa  ,  bsi  =  0.235  nm, 

and  vsi  =  0.218  [5],  a  least  squares  fit  of  Eqn.  1  gives  ~60  dislocations  (see  Fig.  1).  Future 
work  with  in  situ  TEM  will  more  explicitly  give  the  number  of  slip  planes  and  their  back 
stress. 


Figure  1.  BF  TEM  of  a  173  nm  diameter  Si  nanosphere  a)  before  and  b)  after  TEM  in 
situ  compression,  with  the  arrow  in  (b)  highlighting  evidence  of  residual  plasticity,  c) 
Flow  stress  as  a  function  of  nanosphere  diameter. 

Si  nanospheres  and  nanotowers  were  also  fractured  in  situ  with  a  quantitative 
measurement  of  the  fracture  stress.  Indents  were  run  in  displacement  control,  with  10 
nm/s  loading  and  unloading  rates  and  a  5  sec  hold  at  the  peak  load.  Five  spheres  (d=l  13- 
349  nm)  and  three  towers  {d=2'i\-A\5  nm)  showed  signs  of  fracture  at  failure.  Figure  2a,b 
gives  a  representative  example  of  the  fracture  of  a  349  nm  diameter  Si  nanosphere  on  an 
AI2O3  substrate.  To  determine  the  fracture  toughness  of  the  spheres  and  towers,  a  work 
per  unit  fracture  area  method  was  used  since  it  accounts  for  dislocation  plasticity  and 
closely  matches  finite  element  modeling  of  the  fracture  toughness  in  Si  nanotowers  [6]. 
Fig.  2c  summarizes  the  fracture  toughness  for  the  nanospheres  and  nanotowers  along 
with  ex  situ  nanosphere  results  from  Mook  et  al.  [4].  Here,  the  (d  can  be  explained 


theoretically  by  balancing  the  crack  extension  force,  Guc,  with  a  dislocation  resistive 
force,  Fr,  based  on  a  dislocation  line  tension.  With  V5,=0.218  and  taking  yield  strength  as 
the  maximum  value  from  Fig.  Ic,  the  fracture  toughness  becomes 

K,^  =  ^-^/uNbd-^l^  .  (2) 

y 

Eqn.  2  is  shown  in  Fig.  2c  as  a  solid  line,  with  the  same  number  of  dislocations  (A/=60)  as 
determined  in  Fig.  1.  As  with  the  strengthening  argument  proposed  earlier,  this  result 
suggests  that  dislocations  trapped  in  small  volumes  generate  large  back  stresses  and 
increase  the  yield  strength.  In  addition,  the  piled  up  dislocations  will  act  to  shield  the 
crack  tip  and  inhibit  crack  growth,  thus  increasing  the  fracture  toughness. 


Figure  2.  BF  TEM  of  a  349  nm  diameter  Si  nanosphere  a)  before  and  b)  after  TEM  in 
situ  compression,  with  (b)  capturing  a  fracture  event,  c)  Fracture  toughness  as  a  function 
of  nanosphere  diameter. 


A  third  aspect  of  the  deformation  of  Si  in  small  volumes  that  was  pursued  was  the 
role  of  phase  transformation  plasticity.  The  literature  discussing  the  size  effect  on  the  Si 
phase  transformation  process  [7-10]  is  much  more  limited  (and  inconclusive)  than  bulk 
studies.  In  the  present  data,  the  slope  of  the  unloading  curve  for  TEM  in  situ  compression 
of  Si  nanospheres  was  found  to  decrease  abruptly  near  the  end  of  the  unloading,  creating 
an  elbow  (as  shown  in  Fig.  3a).  The  location  of  the  elbow  can  be  defined  as  the  point 
where  the  unloading  curve  deviates  from  the  elastic-plastic  power  law  fit.  This  type  of 
feature  is  commonly  seen  during  the  indentation  of  silicon  and  has  been  attributed  to  the 
transformation  of  the  high  pressure  p-Sn  phase  (Si  II)  of  Si  to  an  amorphous  phase  (a-Si) 
[I  I- 1 3].  Knowing  the  contact  area  between  the  tip  and  the  sphere,  the  pressure  at  which 
this  event  occurs  in  the  212  nm  diameter  nanosphere  is  0.90  GPa.  This  is  less  than  1/6  of 


the  transition  pressure  for  bulk  Si(lll)  [12]  and  suggests  an  inereased  stability  of  Si  11. 
At  even  smaller  volumes,  the  elbow  in  the  unloading  slope  disappears  (see  65  nm 
diameter  sphere  in  Fig.  3b).  This  represents  a  transition  in  the  deformation  mechanism  to 
a  dislocation  plasticity  dominated  response  below  a  critical  length  scale  and  is  currently 
being  investigated  in  more  detail. 
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Figure  3.  Load-displacement  data  for  TEM  in  situ  compression  of  Si  nanospheres  with 
diameters  of  a)  212  nm  and  b)  65  nm.  The  solid  red  line  in  each  is  an  elastic-plastic 
unloading  power  law  fit. 


III.  Si-SiC  core-shell  nanocomposites 

A  second  component  of  this  research  was  in  the  development  of  nanocomposite 
materials  based  on  the  small  scale  strengthening  and  toughening  mechanisms  described 
earlier.  The  HPPD  technique  was  used  to  deposit  nanocrystalline  (<20  nm  grain  size), 
3C-SiC  coatings  with  thicknesses  of  100-200  nm  on  Si(lll)  nanotowers.  The  Si(lll) 
towers  were  vapor-liquid-solid  grown  using  chemical  vapor  deposition  by  collaborators 
at  NIST  [14].  Figure  4  gives  representative  images  of  the  Si  towers  before  and  after 
deposition,  with  the  towers  attached  to  their  original  growth  substrate  in  each  case.  Based 
on  the  DF  TEM  imaging,  the  composite  structures  shows  the  conformal  nature  of  the 
coating,  with  good  evidence  that  the  appropriate  adhesion  has  been  achieved. 


Figure  4.  Examples  of  Si(l  1 1)  nanotowers  from  batch  B  a)  before  (BF  TEM)  and  b)  after 
deposition  of  the  SiC  coating  (DF  TEM).  The  parallel  lines  rurming  down  the  length  of 
the  Si  in  (b)  are  from  thickness  fringes. 

The  composite  towers  were  then  heat  treated  using  rapid  thermal  annealing  (RTA) 
for  20  sec  at  900-1200°C  in  an  Ar  environment.  Due  to  the  mismatch  in  the  coefficients 
of  thermal  expansion  for  the  Si  core  and  SiC  shell,  the  Si  core  was  left  in  residual 
compression  following  the  RTA  treatment.  This  compressive  stress  was  measured  using 
DF  TEM  strain  contours,  confocal  Raman  microscopy  (CRM),  and  atomic  force 
microscopy  (AFM)  as  shown  in  Figures  5  and  6.  Oscillations  in  DF  TEM  contrast  (see 
Fig.  5a)  originate  from  the  base  of  the  tower  and  show  axial  symmetry.  This  is  consistent 
with  an  elastic  strain  profile  for  a  capped  annular  coating,  where  the  highest  strain  is  at 
the  top  of  the  Si  tower  and  relaxes  at  the  base  due  to  less  confinement.  CRM 
measurements  used  a  5 14.5  nm  Ar  laser  focused  perpendicularly  onto  the  (1 1 1)  surface  of 
the  Si  nanotower.  Measurements  taken  after  RTA  heat  treatments  showed  a  significant 
blue  shift,  with  the  strongest  shift  for  the  1200°C  treatment  (see  Fig.  5b).  Assuming  the 
Si  towers  are  under  a  biaxial  stress  from  the  CTE  mismatch  with  the  SiC  coating,  this 
blue  shift  represents  a  compressive  stress  of  0.8-2. 8  GPa  (see  Fig.  5c).  While  the 
magnitude  of  this  stress  is  3-5  times  larger  than  an  elastic  estimate  for  concentric 
cylinders,  the  trend  of  increasing  stress  with  increasing  RTA  temperature  was  verified  by 
cross  sectioning  the  towers  along  the  radial  axis  and  measuring  the  stress  relaxation  of  the 
Si  core.  Figure  6a, b  shows  the  composite  surface  that  was  cross  sectioned  with  the 
focused  ion  beam.  This  produces  a  nominally  flat  surface,  with  the  contrast  in  Fig.  6b  due 
to  the  compositional  difference  between  the  core  (Si)  and  shell  (SiC).  Following  RTA 
treatment,  the  height  of  the  extruded  Si  core  was  measured  using  AFM  contact  mode 


imaging  (Fig.  6c, d).  With  increasing  RTA  treatment  temperature,  there  is  an  increased 
volume  of  Si  extruded  due  to  the  higher  CTE  stresses  on  the  core. 
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Figure  5.  Evidenee  of  residual  stress  in  Si  eore,  showing  a)  DF  TEM  strain  eontours,  b) 
blue  shift  in  the  Raman  Si  response  (normalized,  Lorentzian  fittings)  for  RTA  treated 
(900-1200°C)  composite  towers,  and  c)  the  biaxial  stress  of  the  Si  core  as  a  function  of 
RTA  treatment  temperature. 


Figure  6.  Additional  evidence  of  residual  stress  on  FIB  cut  composite  Si  towers,  showing 
a,b)  the  FIB  cut  surface  exposing  the  eore-shell  structure  (two  different  magnifications), 
c)  AFM  contact  mode  imaging  of  a  1200°C  RTA  treated  tower  with  an  extruded  Si  core, 
and  d)  extruded  core  height  as  a  funetion  of  RTA  treatment  temperature. 


Considering  the  inereased  stability  of  Si  II  in  small  volumes  (as  shown  in  Seetion 
II)  together  with  the  high  eompressive  stresses  found  in  the  Si-SiC  eore-shell  composites 
(Section  III),  a  new  toughening  mechanism  for  Si-SiC  nanocomposites  was  proposed  as 
follows.  First,  Si  I  nanospheres  with  diameters  in  the  range  of  100-1000  nm  are  dispersed 
in  a  SiC  matrix  and  heat  treated  to  1 100-1200°C.  Based  on  Fig.  5c,  this  will  leave  the  Si 
nanospheres  under  a  compressive  stress  of  close  to  3  GPa.  The  composite  is  then  loaded 
to  a  stress  exceeding  the  critical  stress  needed  for  the  Si  I^II  transition  (~8-10  GPa). 
This  can  be  done  either  under  service  conditions  or  by  pre-stressing  the  composite.  Upon 
unloading,  the  nanospheres  will  remain  as  Si  II  since  the  residual  compressive  stress  is 
higher  than  the  Si  II— >I  transition  pressure.  When  a  crack  forms  in  the  composite,  the 
tensile  stress  at  the  crack  tip  will  destabilize  the  Si  II  and  cause  a  transformation  back  to 
Si  I.  The  volume  expansion  during  this  transformation  will  then  act  as  a  compressive 
stress  to  close  the  approaching  crack,  in  same  way  transformations  in  partially  stabilized 
zirconia  (PSZ)  have  been  shown  to  cause  crack  closure  [15,  16].  If  the  tensile  stress  at  the 
crack  tip  does  not  transform  the  Si  II,  the  Si  II  nanospheres  will  act  as  either  ductile 
inclusions  for  crack  pinning  or  form  voids  due  to  the  Si  I— >II  contraction  that  will  serve 
as  crack  arrest  points.  Thus,  by  controlling  the  length  scale  of  the  Si  nanospheres,  Si-SiC 
nanocomposites  can  be  toughened  by  a  combination  of  phase  transformation  and  ductile 
phase  reinforcement  mechanisms. 

IV.  MD  simulations 

Molecular  dynamics  simulations  were  used  to  investigate  the  mechanical 
properties  resulting  from  the  compression  of  silicon  nanospheres,  with  a  particular 
emphasis  on  determining  the  cause  of  the  experimentally  observed  hardening  behavior. 
The  yielding  properties  observed  within  these  modeled  silicon  nanospheres  showed  a 
high  dependence  on  changes  in  the  atomic  potential,  temperature,  orientation  and  sphere 
size.  These  variations  allowed  for  the  study  of  numerous  yielding  behaviors,  including 
phase  transformations,  dislocation  nucleation  and  motion,  and  stacking  fault  formation 
and  growth. 


Figure  7.  (a)  Cross  sectional  image  of  a  10  nm  sphere  at  0  K  using  the  Tersoff  potential 
showing  that  the  center  has  transformed  into  P-Sn  (yellow),  (b)  Compression  of  a  34  nm 
diameter  sphere  at  300  K  using  the  Stillinger- Weber  potential  reveals  multiple  dislocation 
loops  that  are  seen  to  nucleate  from  the  contact  areas. 


By  relating  the  atomistic  yielding  mechanisms  to  the  measured  hardening 
behavior,  it  was  found  that  the  hardening  behavior  was  independent  of  the  pressure 
induced  transformation  to  the  P-Sn  phase.  Strong  indications  were  also  found  that  this 
hardening  could  be  the  result  of  dislocation  interactions.  To  better  investigate  this,  a 
million  atom  simulation  was  performed  to  maximize  the  number  of  dislocations.  This 
revealed  numerous  dislocations  that  were  able  to  interact  resulting  in  a  short  period  of 
hardening.  The  main  limitation  observed  was  that  the  dislocations  quickly  reach  the 
sphere's  surface  and  disappear  due  to  the  sphere's  lack  of  an  oxide  layer.  Currently,  work 
is  being  done  looking  at  how  the  presence  of  this  oxide  would  affect  the  dislocation 
behavior. 


V.  Conclusions 

To  date,  the  experiments  performed  have  made  substantial  progress  towards 
achieving  the  objectives  of  this  project.  Inverse  length  scale  dependent  relationships  for 
strength  and  toughness  in  Si  were  identified  based  on  dislocation  pile-up  and  crack  tip 
shielding  mechanisms,  respectively.  In  addition,  a  transition  was  identified  in  the 
deformation  of  Si  under  anisotropic  loading  below  a  critical  size,  showing  increased 
stability  of  Si  II  and  the  transition  at  a  critical  length  scale  to  dislocation  dominated 
plasticity.  By  imposing  a  highly  compressive  stress  on  the  Si  core  in  Si-SiC  composites. 


the  phase  stability  of  Si  II  was  applied  to  a  new  toughening  mechanism  in  Si-SiC 
composites.  Finally,  MD  simulations  verified  the  presence  of  dislocation  based  plasticity 
during  the  deformation  of  small  volumes  of  Si.  Overall,  these  results  demonstrate  the 
importance  of  nanoscale  confinement  and  localized  stress  in  the  design  of  mechanically 
robust  nanocomposites. 

Currently,  experiments  are  underway  to  confirm  the  toughening  mechanism 
proposed  in  Section  III  through  a  collaboration  with  EMPA  Swiss  Federal  Laboratories. 
These  experiments  involve  the  SEM  in  situ  compression  of  RTA  treated  Si-SiC 
composite  towers  to  a  pressure  above  the  Si  II  transition.  Upon  release,  we  expect  to  see 
the  Si  II  phase  locked  into  place  and  will  be  able  to  verify  this  by  the  absence  of  an  elbow 
in  the  load-displacement  unloading  curve  and  CRM  analysis.  In  addition,  we  plan  to  FIB 
cut  precracks  in  the  composite  towers  and  compress  the  towers  in  situ  to  observe  the 
crack  interaction  with  the  compressively  stressed  Si  core.  Finally,  we  are  beginning  HR 
TEM  experiments  with  in  situ  deformed  Si  nanotowers  in  the  labs  of  Prof  Ikuhara  to 
identify  the  types  and  arrangements  of  dislocations. 
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